Parametric form of QCD travelling waves 
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We derive parametric travelling-wave solutions of non-linear QCD equations. They describe the 
evolution towards saturation in the geometric scaling region. The method, based on an expansion 
in the inverse of the wave velocity, leads to a solvable hierarchy of differential equations. A universal 
parametric form of travelling waves emerges from the first two orders of the expansion. 



1. Since the last years, progress is being made in the understanding of the QCD evolution in rapidity towards 
saturation, i.e. the high-density limit of weekly interacting partons probed at high energy. In particular, a link has 
been established Q between travelling-wave solutions of non-linear equations and geometric scaling of the gluon 
distribution function N(L, V), where L = logfc 2 is the transverse momentum variable and Y the rapidity. 

In the mean-field approximation we will consider here, J\f obeys the nonlinear Balitsky-Kovchegov (BK) evolution 
equation 

"lh: d Y N = x {-d L )N -N 2 , (1) 

o ■ 

where, for convenience, rapidity Y is measured in units of a " Nc the strong coupling constant, kept fixed. 

(N : 



a : 



x(7) = 2<Kl)-<K 7 )-V(l-7) (2) 

is the Balitsky Fadin Kuraev Lipatov (BFKL) kernel Q . 

The identification of Eq. JJJ , in the diffusive approximation of the kernel, with the Fisher and Kolmogorov-Petrovsky- 

Piscounov (F-KPP) non-linear equation || infers Q the existence of asymptotic travelling wave solutions in rapidity, 

£Nj , see Fig^ The mechanism for the formation of travelling waves is more general and applies also Q] to the BK 

equation Eq.QJ. Starting with the initial condition Af(L, Yn) ~ exp (—jqL), travelling wave solutions Af(L—c Y) are 

' ' , formed during the F-evolution with the following features [3 : 
^ , 

C"*** ■ • "pushed front": 70 < j c ; c= x{l)/li with 7 = 70 

m ■ 

' • "pulled front": 70 > 7c! c = x(l)/'7j with 7 = 7 C , 

in ! 

where, for the BK equation, 7 C = .6275 is the critical anomalous dimension, solution of x(7c)/7c = x'ilc) — c, where 
c = 4.883 is the critical velocity and 70 = 1 > j c (color transparency). For the F-KPP equation, the corresponding 
values are j c = 1 and c = 2. These powerful results come from a "locking" of the solutions of the linear part of the 
kernel due to the non-linearities. For the "pulled fronts", it has universal properties, i.e. the locking is independent 
Q-i! of the form of non linearities and initial conditions. For the "pushed fronts" it is not so. 

However, there are limitations of these solutions. They may require very large ranges in Y and L to appear. 
Moreover, while the wave structure is clear, the form of the front deduced from the linear kernel does not extend to 
1-^ ■ the region where non-linear effects are important. In fact, only very few information on the non- linear term has been 
^ , used. The goal of the present paper is to present a different method allowing to analytically derive travelling wave 
solutions taking into account the non-linear terms of the equation. 

The values of 7 and c coming from the knowledge of asymptotic properties are considered as input parameters for 
travelling wave solutions of the full non- linear equation, and one uses the fact, already noticed in the literature for 
the F-KPP equation @], that l/c can be taken as a small expansion parameter of the full solution. 

2. Let us start with Eq.JU, where the integro-differential operator \ will be defined with the help of a limited 

expansion around the value of 7 relevant for the travelling wave solutions, i.e., either 70 or 7 C depending whether we 
are in a "pushed" or "pulled" front situation. 

X (-9l) = X(7)l + x'm-dL - 7l) + ix"(l)(-d L - 7I) 2 + iX i3) (l)(-d L - 7 1) 3 + ■ ■ • • (3) 

We shall consider truncations of in the number of derivatives to some value P. We are looking for a solution in the 
vicinity of the selected anomalous dimension 7. For instance the diffusive approximation equivalent to the F-KPP 
equation is obtained with P = 2. The validity of the approximation will be discussed. 
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A generic travelling wave (i.e. geometric scaling) solution has the form M{z = bL~ aY), with a, b fixed parameters. 
Inserting it into Eq.Q), the partial differential BK equation Q turns into an ordinary differential equation 



A Q N - + (a - A lb )^- + E(-1) P ^ 6P S = 

p=2 ^ ' 



where 
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With suitable redefinitions 



a 



z/c; 



one obtains the equation 
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with A = \J A$A 2 . 

Eq. JJJ provides an expansion of the nonlinear differential equation in terms of 1 /cP terms with the 1 /c term missing. 
Limiting J2J to the three first terms (A — > 0) gives the equation that one obtains [2( for travelling waves in the F-KPP 
equation 



dtu(t, x) = d 2 u(t, x) + u(t, x)(l — u(t, x)) 



(8) 



by inserting the form u(t, x) = U(s — x/c—t). Hence in the expansion around the F-KPP equation, the critical velocity 
gets renormalized from its original value to c ~ 2 by the redefinitions ©. We see that adding higher derivatives around 
7 up to rank P leads to a hierarchy of terms in l/c p , p > 3 with c > 2. 

Taking 1 / c as the small parameter of an expansion we are looking for an iterative solution 



h{s) = flO + ^h 2 + E ~r 



p>3 



K = -- u(s). 



Inserting it into Eq.lJJJ) translates into a hierarchy of equations 

h' + hi - 1/4 
h' 2 + 2h h 2 + h' ' 



tin + 2h a h 3 
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h' A + 2h h 4 + hj + h' 2 ' - \A4,ti Q "' = 



h' + 2h h 5 + 2h 2 h 3 



\A 3 ti 2 " 



(9) 



where we have restricted (for sake of simplicity) the hierarchy of equations to 0(l/c 6 ). An iterative solution can easily 
be found with initial conditions being appropriately chosen, i.e. ft,o(±oo) = ±i and hi^o(±oo) = hi(0) = for higher 
orders. Note that if h(s) is solution, also h(s + s ) is solution. 

The system is iteratively but fully solvable. One first solve the only non-linear equation of the hierarchy, obtaining 
ho = fth(§). Using the property 



ds 



h n (s) + 2h h n {s) 



1 



ch 2 {s/2) ds 



[ch 2 (s/2) h n (s)] , 



(10) 



all the other linear equations reduce to simple integrations of functions defined recursively from the hierarchy. 
The solution of the first three terms of the expansion is 
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(11) 




FIG. 1: Universal travelling wave. The function u(x,t) = U(s — x/c — t) is represented in the critical case (c=2) for three 
different times. Continuous lines: U(x/2 — t); Dashed lines: First C(l/c°) term only. The wave front connecting the regions 
u = 1 and u — travels from the left to the right as t increases. That illustrates how, in terms of QCD properties, the 
"saturation" region invades the "transparency" region. Note that the first order correction is already small. 



This can be compared with the hierarchy of equations obtained for the F-KPP equation JSJ) 
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2h h 2 
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where the expansion parameter is now 1/c 2 . The solution of the first three equations 8] of l|12l) is 



(12) 
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l + e s 
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(e s + l) 2 



(13) 

The comparison of the sets of equations 10 and (|12H and their solutions (|llfl and (| 1 31) calls for comments. 

i) The two first terms of the expansion in the small parameter of order 1 and 1 /c 2 are the same, and then give 
an "universal parametric form" to the travelling wave solution, see FigQ The form is unique while the parameters 
values may change by a redefinition © of the coefficients Ai. 

ii) The non- universal dependence on the kernel appears at order 1/c 3 . It can stay small provided the coefficient A% 
and the followers remain bounded. This implies to keep the expansion of the kernel limited to the region around the 
chosen value of 7 and in particular in the vicinity of 7 C in the critical case c = 2. 

iii) The truncation of the expansion of the kernel is required for the BK equation equipped with the BFKL kernel. 
Indeed, if considering the full analytic expansion to define the coefficients Ai (i.e. P — ► 00), they are infinite since 
they correspond to the expansion of the BFKL kernel near 7 = where it possesses a pole singularity. We interpret 
this phenomenon as a breaking of the travelling wave property in the deep saturation region of the BK equation. In 
fact, both in the deep saturation region (u ~ 1 in Fig.JU) and in the leading edge of the wave (it ~ 0) we expect the 
exact scaling properties not to be valid Scaling is expected in between, that is in the "interior" of the wave. 

3. Using the redefinitions © , and limiting ourselves to the universal part, we obtain 
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where 



Q 2 s (Y)=exp(n 2 Y) 



plays the role of the saturation scale. In (|14ll 
H14JI only depends on the function 



we have defined y,\ 
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The solution 
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(16) 



where one writes W 2 = e Y , Wo = k^^ 1 ^ 2 . Hence e s has the typical structure of a forward Regge amplitude. While the 
first form exhibits the geometric scaling structure, the second one suggest an interpretation of (|14|) as an unitarization 
of the amplitude (fTTT^l . 

Hence, we have found the possibility of assessing a specific parametric form to the geometric scaling solutions of 
the BK equation, through a redefinition of the parameters fii,^. The form l|14l) . with the possibility of extending the 
parametrization to higher orders by a solution of the hierarchy may help establishing a link between the theory 
and the phenomenological observation of geometric scaling Q. This is the subject of a subsequent work 

On a more theoretical level, it seems feasible to extend the method to QCD evolution equations beyond the "mean 
field" BK equation Q). In particular the extension to the stochastic versions of QCD evolution equations could 
be treated |9( by averaging the parametric form over the distribution of velocities c generated by the noise. 

In conclusion, the fruitful properties of travelling- wave solutions of non-linear evolution equations, both through 
the asymptotic limit and through their expansion in terms of the inverse velocity give significant tools to deepen our 
understanding of QCD in the high density limit and emphasize once more the interesting connections with statistical 
physics. 
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